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Abstract 

Waveguiding  in  a thin-film  helicoidal  bianistropic  medium  (TFHBM)  layer  is  investigated. 
A dielectric  TFHBM  layer  bounded  by  isotropic  dielectric  half-spaces  is  shown  to  support 
guided  wave  propagation  with  guide  wavenumbers  dependent  on  the  direction  of  signal  prop- 
agation, thus  signalling  potential  use  as  a space-guide.  The  modal  fields  and  power  transmis- 
sion distributions  associated  with  the  guided  modes  in  the  proposed  TFHBM  interconnects 
are  detailed. 


1.  Introduction 

Implementation  of  optoelectronic  devices  requires  the  development  of  optical  interconnects 
which,  in  addition  to  providing  effective  signal  transmission,  must  be  simple  to  fabricate  on 
integrated  circuitry.  In  this  paper,  we  present  a theoretical  study  which  indicates  that  dielectric 
thin-film  helicoidal  bianisotropic  mediums  (TFHBMs)  are  very  suitable  for  realizing  optical  in- 
terconnects. In  fact,  the  adoption  of  dielectric  TFHBM  interconnects  may  result  in  efficient  use 
of  semiconductor  real-estate  in  electronic  chips. 


2.  Theory  in  Brief 

Suppose  a linear,  dielectric  TFHBM  layer  completely  fills  the  region  \z\  < D/2,  while  the 
halfspaces  z < —D/2  and  z > D/2  are  filled  by  an  isotropic  dielectric  medium  whose  rela- 
tive permittivity  scalar  at  angular  frequency  w is  denoted  by  er(u;).  The  constitutive  relation 
D(r,w)  = e0g(z,u>)  • E(r,u;)  of  the  TFHBM  layer  contains  the  relative  permittivity  dyadic 

g(z, u>)  =gz{z-%)  * 5,(x)  * £e/(u>)  • §F(x)  • g(z  - f ) . (1) 

Here,  the  reference  relative  permittivity  dyadic 


£ef(w)  = fa(«)n*n*  + e&(w)u*u*  + ec(cj)Uj,Uj, , 


(2) 


where  ux>y>z  are  the  cartesian  unit  vectors  and  eQ,t,c(w)  are  frequency-dependent  scalars.  The  ro- 
tational non-homogeneity  (along  the  z axis)  of  a structurally  right-handed  TFHBM  is  expressed 
by  the  dyadic 


( 7TZ\  { 7TZ\ 

§^{z)  = (u*ux  + UyUy)  cos  f — J + (uyux  - uxuy)  sin  ( — 1 + uzuz 


(3) 
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with  2fJ  as  the  structural  period.  The  so-called  angle  of  rise  x appears  in  the  tilt  dyadic 


iyix)  = («xUx  + Uzuz)  cos  x + (t^Ux  - uxu*)  sinx  + Uj,Uj, ; (4) 

typically  x 20°  for  the  sculptured  thin  films.  Guided  wave  propagation  is  ensured  if  er  < 
min{ea,  ej,,  ec},  with  both  mediums  assumed  non-dissipative  at  the  frequency  of  interest. 

Knowing  the  constitutive  relations  of  the  chosen  TFHBM  layer,  we  can  determine  the  guide 
wavenumbers  which  enable  guided  wave  propagation.  A specific  guided  wave  mode  can  be 
delineated  with  the  following  equations: 


E(r)  = exp  [in  (x  cos  ip  + ysini/>)]  e(z,n,ip)  1 

H(r)  = exp  [*k  (x  cos  ip  + y sin  ip)]  h (z,n,ip)  ( ’ 


— oo  < z < oo . 


(5) 


Here,  the  angle  ip  denotes  the  propagation  direction  = ux  cos  ip  + uysm  ip,  while  k is  the 
modal  guide  wavenumber  whose  values  have  to  be  determined. 

The  leakage  fields  accompanying  a guided  wave  mode  are  represented  by 

E(r)  = ( bss  + 6pp_) exp  [ik_  • (r  + § uz)]  1 < _£) 

H(r)  = |(6ap_  -6ps)exp  [ik_  • (r  + fuz)]  J ’ 2 

in  the  lower  halfspace,  and 

E(r)  = (cas  + cpp+)exp  [ik+  • (r  - § u*)]  1 D 

H(r)  = |(c*P+  -Cps)exp  [zk+  • (r-  §uz)j  J ’ Z ~ 2 ’ ^ 

in  the  upper  halfspace,  with  bs , bp , cs , and  Cp  as  the  amplitudes  of  the  perpendicular-  and  the 
parallel-polarized  components.  The  various  vectors  introduced  in  (6)  and  (7)  are  given  by 

s = u2xu<]  P±  = T [l  - (/c/fc)2]  U/+  ( n/k)uz  , k±  = ku£  ± (k2  - k2)  1/2  uz  , (8) 

where  k = /cov/e7,  V = Vo/s/^,  K = aye 0/i0  = 2n/X0  is  the  free-space  wavenumber,  A0  is  the 
free-space  wavelength,  and  rj0  = \Z(J.0/e0  is  the  intrinsic  impedance  of  free  space.  Guided  wave 
propagation  is  possible  only  if  k > k;  otherwise,  energy  launched  into  the  TFHBM  layer  must 
leak  into  the  two  halfspaces. 

On  substituting  the  constitutive  relations  and  the  field  expressions  above  into  the  time- 
harmonic  Maxwell  curl  equations  V x E(r)  = zwB(r)  and  V x H(r)  = —iui€0e(z,uj)  • E(r),  a 
4x4  matrix  ordinary  differential  equation  emerges  for  the  TFHBM  layer.  Its  solution  requires 
the  prescription  of  boundary  values  through  (6)  and  (7).  Our  interest  lies  in  determining  pairs  of 
( k , ip)  such  that  not  all  of  the  coefficients  6S,  6P,  c9  and  Cp  are  null-valued;  thereby  the  dispersion 
equation  is  obtained.  For  guided  wave  propagation,  values  of  k denoted  by  «£,,  (r  = 1,  2,  3,...), 
that  satisfy  the  dispersion  equation  have  to  be  numerically  determined,  the  roots  being  indexed 
by  the  integer  r in  descending  order  of  their  magnitudes. 


3.  Guide  Wavenumbers 

We  implemented  the  foregoing  procedure  using  the  C programming  language  and  the  IMSL 
C numerical  library  subroutines  for  complex  linear  algebra.  The  wavelength  A0  was  fixed  at 
600  nm  for  all  calculations.  We  tested  our  computer  program  for  the  case  of  an  isotropic, 
homogeneous,  dielectric  slab  waveguide.  Analytical  solutions  to  the  dispersion  equation  of  this 
simple  waveguide  are  well-documented  - see,  e.g.,  [1].  Setting  er  = 1,  we  simulated  homogeneity 
and  isotropy  by  choosing  eQ  = e6  = ec  and  taking  the  limit  I/O,  — ► 0.  The  roots  Kr  that  we 
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obtained  corresponded  exactly  with  the  analytical  results.  Furthermore,  the  calculated  mode 
shapes  and  power  transmission  characteristics  of  the  waveguide  also  matched  the  expected  power 
transmission  and  mode  shape  plots. 

Now  let  us  proceed  to  the  proposed  TFHBM  interconnect.  Most  calculations  were  made  with 
{e0  = 3.8,  ej,  = 4.6,  ec  — 3.0},  in  accordance  with  data  from  [2],  In  general,  the  guide  wavenum- 
bers show  a strong  dependence  on  the  propagation  direction  (see  [3,4]  for  more  details).  This  is 
illustrated  in  Figure  1,  where  the  guide  wavenumbers  are  indicated  for  various  ip  for  a specific 
TFHBM  interconnect.  Thus,  the  proposed  TFHBM  interconnect  functions  as  a space-guide 
through  which  signals  can  be  simultaneously  transported  in  different  directions  with  different 
phase  velocities.  This  feature  emerges  from  the  anisotropic  and  non-homogeneous  nature  of 
TFHBMs,  and  may  be  exploited  for  efficient  use  of  semiconductor  real-estate  in  optoelectronic 
circuitry. 
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Figure  1:  Roots  of  dispersion  equation  for  directions  of  propagation  denoted  by  the  angle  ip-, 
A0  = 600  nm,  ea  = 3.0,  ej,  = 4.6,  ec  = 3.8,  x = 30°,  D — 200  nm,  D = 8 0 = 1600  nm. 


Independently  of  all  parameters,  the  guide  wavenumbers  kL  for  propagation  directions 
and  — Uf  are  the  same.  When  the  TFHBM  layer  consists  of  an  integral  number  of  periods  (i.e., 
the  ratio  D /fl  is  an  even  integer),  the  additional  relation  holds.  This  arises  because 

all  three  principal  axes  of  g (z,  to)  rotate  through  an  integral  number  of  turns  between  the  planes 
z = —D/2  and  z = D/2,  thereby  imposing  a symmetry  constraint. 

The  variability  of  k ^ with  ip  is  most  pronounced  around  the  lower  values  of  kt^  (where 
solutions  of  the  dispersion  equation  are  more  widely  spaced).  Additionally,  the  directional- 
dependence  of  persists  for  smaller  values  of  D/D,  including  D/D  « 1.  Parenthetically,  we 
also  studied  the  directional-dependence  of  the  guide  wavenumbers  for  a locally  uniaxial  TFHBM 
layer  with  {e0  = ec  = 3.8,  ej,  = 4.6};  The  guide  wavenumbers  kL  still  exhibit  a dependence  on  ip, 
but  the  dependence  is  weaker  than  for  the  biaxial  case  illustrated  in  Figure  1. 

The  number  density  of  guided  wave  modes  is  less  when  D is  small.  The  mode  number  density 
appears  to  be  predominantly  determined  by  the  overall  thickness  D,  and  is  largely  unaffected 
by  the  half-period  D.  Thus,  the  availability  of  guided  wave  modes  can  be  tailored  by  properly 
choosing  the  layer  thickness  D. 

In  TFHBM  interconnects,  has  an  upper  bound  which  varies  with  ip.  For  instance,  kL  < 
2.049  k0  for  ip  — 0°,  whereas  < 1.981  k0  for  ip  = 90°,  in  Figure  1.  The  upper  bound  decreases 
monotonically  as  ip  increases  from  0°  to  90°.  The  upper  bound  on  varies  with  x,  ip,  D,  and 
D for  given  eQ,  e*,,  and  ec. 


184 


4.  Space— Guide  Modes 

The  modal  fields  and  power  transmission  associated  with  the  guided  wave  modes  were  also 
studied.  The  time-averaged  power  flow  in  the  propagation  direction  is  given  by  Pe{z)  = 
luf  • R e[e(z,n,ip)  x h.*(z,K,4>)],  where  the  asterisk  denotes  the  complex  conjugate.  Detailed 
numerical  study  of  the  modal  fields  and  power  distributions  revealed  that  each  mode  of  propaga- 
tion inside  a TFHBM  interconnect  can  be  classified  into  one  of  two  groups:  hybrid  electric  (HE) 
and  hybrid  magnetic  (HM).  The  modes  are  hybrid , because  electric  and  magnetic  field  compo- 
nents are  present  in  all  directions,  along  the  axial  (uz)  direction,  as  well  as  in  the  longitudinal 
and  perpendicular  directions  in  the  xy  plane.  This  is  unlike  the  modes  in  an  isotropic,  dielectric, 
planar  interconnect,  wherein  the  modes  are  either  transverse  electric  (TE)  or  transverse  mag- 
netic (TM).  Another  distinction  between  the  HE  and  HM  modes  in  the  TFHBM  interconnect 
and  TE  and  TM  modes  in  the  isotropic  interconnect  is  that  there  appears  to  be  no  apparent 
ordering  to  the  occurence  of  the  HE  and  HM  modes,  while  the  TE  and  TM  modes  alternate 
with  r = 1, 2, 3, ... . 

For  both  the  HE  and  the  HM  modes,  the  power  transmission  distributions  Pe(z)  are  quite 
similar  to  those  of  the  TE  and  the  TM  modes,  respectively.  In  fact,  for  propagation  in  any 
direction,  it  is  possible  to  order  the  guided  modes  HEn  and  HMn,  (n  = 1,2,3,...),  based  upon 
the  similarity  of  Pt{z),  respectively,  to  Pe{z)  for  TEn  and  TMn  modes.  Pi{z)  for  a given  mode 
(HEn  or  HMn)  does  not  vary  much  with  respect  to  1/7  in  the  space-guide.  Regarding  modal  field 
plots,  however,  there  are  distinct  differences  between  the  HE  and  the  HM  modes. 

The  variation  of  ej.  = e • u e with  respect  to  z for  the  HEn  mode  is  similar  in  all  directions; 
and,  in  general,  the  ej_  vs.  2 curves  for  all  HEn  modes  resemble  those  for  a TE„  mode  in 
an  isotropic  interconnect.  Thus,  all  HEn  modes  propagating  in  any  direction  in  a space- 
guide  have  similar  modal  characteristics.  However,  the  h±  vs.  2 plots  for  the  HMn  modes  do 
not  display  these  characteristics.  Not  only  are  the  h±  vs.  z plots  for  an  HMn  mode  different 
from  that  of  the  TMn  mode,  but  also  the  h±  vs.  z plots  are  distinctly  dissimilar  for  the  various 
propagation  directions.  Thus,  the  dielectric  anisotropy  and  non-homogeneity  of  the  space-guide 
impart  more  significant  directional  dependence  to  the  mode  shapes  of  the  HM  modes  and  less 
to  the  HE  modes. 

Clearly,  the  HEn  mode  launched  in  one  direction  will  not  interfere  with  the  HEn  mode 
launched  in  some  other  direction;  and  the  same  holds  true  for  any  HMn  mode,  at  least  for 
small  values  of  n.  Indeed,  several  HE  and  HM  modes  of  low  order  can  be  launched  in  different 
directions,  while  taking  care  that  their  guide  wavenumbers  are  all  different.  The  space-guide 
concept  is  thus  well-founded.  Obviously,  however,  hardware  requirements  will  put  a limit  on 
the  number  of  channels  a TFHBM  space-guide  can  realistically  support  in  actual  circuitry. 
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